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Abstract: This paper deals with multiple -input 
multiple-output orthogonal frequency division 
multiplexing (MIMO OFDM) system. The 
performance of a linear minimum mean square error 
(LMMSE) channel estimator is estimated using power 
delay profile (PDP) technique, hi estimating the 
power delay profile (PDP), all the pilot symbols must 
be used at transmit antennas. The insufficient 
number of samples and null subcarriers causes 
distortions which leads to inappropriate PDP 
estimation. So to obtain the accurate PDP the 
implemented technique holds good. The simulated 
result using the implemented PDP techniques shows 
that LMMSE channel estimation perfor 
approaches the wiener filtering. 



INTRODUCTION: 

Multiple input multiple output is the dominant air 
interface for 4G and 5G broadband wireless 
communications. MIMO combines multiple input 
multiple output which multiplies capacity by 
transmitting different signals over multiple 
antennas. The idea behind OFDM is to divide the 
available spectrum into large number of closely 
spaced subchannels. The combination of MIMO 
and OFDM technology achieves high spectral 
efficiency. By using multiple antennas and 
preceding the data different data streams could be 
sent over different paths. MIMO at higher speeds 
would be most manageable using OFDM 
modulation, because OFDM converts a high speed 
data channel into a number of parallel, lower speed 
channels. The advantage of MIMO OFDM is 
obtaining a spatial diversity which provides 
considerable gain compared to broadband single 
antenna. In MIMO OFDM[l] systems channel 
estimation plays an important role, so it must be 
known accurately to provide considerable gain at 
the receiver. 

When the receiver knows the channel 
statistics, then the mean square error is optimized 
by using LMMSE technique for pilot aided channel 
estimation. The power delay profile estimation is 
implemented to obtain the channel statistics [2] at 
the receiver. OFDM has the advantage of cyclic 
prefix [10] (CP), because of this advantage 
maximum likelihood(ML) estimation schemes are 



used. To obtain the accurate PDP [9] the ML PDP 
estimator requires complex computation. 

By estimating the second -order channel 
statistics i.e, mean delay and root-mean- 
square(RMS) delay spread can improve the 
LMMSE channel estimation performance with 
approximated PDP. Using pilots the channel delay 
parameters are estimated with low complex 
computations .For practical applications like 
WiMax systems the LMMSE estimator [6] with 
PDP is appropriate. The correlation mismatch and 
the error in delay parameters estimation results in 
performance degradation. 

To improve the performance, mismatch 
must be reduced. The PDP technique implemented 
in this paper reduces the mismatch. This 
implemented technique has less computational 
complexity in estimating the PDP, why because 
only the pilots from all transmit antenna ports is 
used. It also reduces the distortion effects caused 
by insufficient number of estimated channel 
impulse response(CIR) samples. The simulation 
results obtained using the implemented technique 
shows the performance of LMMSE channel 
estimations [5] approaches the Wiener filtering. 
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Fig 1: OFDM Transmitter 
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Fig 2: MEMO -OFDM Receiver 

SYSTEM MODEL: 

The system has P transmit and Q receive antenna 
MIMO OFDM system is considered with K total 
subcarriers fe d data subcarriers and K-k d pilot 
subcarriers. The analog signal is converted to 
digital signal and is encoded using the encoder the 
pilots are inserted to output of encoder. Then the 
pilot aided OFDM signal is modulated using QAM 
modulation technique then CP is added to the 
modulated signal to eliminate the interference. The 
length of CP is must be longer than channel delay 
L ch , then only adding CP is beneficial. Adding CP 
is nothing but the performance of IFFT, after IFFT 
the obtained output is the OFDM signal. Now it is 
transmitted to the wireless channel. At the MIMO 
receiver the reverse operation is performed. CP is 
removed FFT operation is performed. The 
corresponding transmit and receive antenna having 
CIRs has same PDP. At the qth receive antenna the 
pilot symbol is represented as 

y q [n p ] = diag(X p) F p h pq + n q (1) 

Where, 

h p , q = [h p , q [n p ,0] ...h pq [n p ,L ch ],0..,0] T 
Which is an L g X\ CIR vector at pth ,qth transmit 
and receive antenna. 

X p = [C p [i 1 ,n p ],C p [i 2 ,n p ],..,C p [i kp ,n p ]] T 
is a pilot vector at the rtpth OFDM symbol 
fort k £ F p andk=l,2..,tf p . 

diag{X p ) is K p X K p diagonal matrix whose entries 
are the K p elements of the vector X p . 
F p is a K p XL g matrix with the (i k ,l)th entry 
[F p ] ik ,l = l/Vtfexp{-j27ri k J/tf} where, i k £ F p 
and 1=0,1, ....,L g - 1. 

n q = complex additive white Gaussian noise vector 
(AWGN) which is having zero mean and 
variance o% at each entry at the q th receiver 
antenna. 



IMPLEMENTED METHOD: 

Using the regularized least squares (RLS) channel 
estimation with a fixed length l g CIR at the (p,q) th 
antenna port can be estimated as 

Ji^^F^ehj'F^diagiX^y^] 
= W RLSp y q [n p ], (2) 

Where, 

£= 0.001 is a small regularization parameter, and 
I Lg is an identity matrix with size L g XL g . 
Due to spars ity of pilot tones in frequency domain 
and the presence of virtual subcarriers makes F p F p 
in (2) ill-conditioned. The ensemble average gives 
the PDP from the estimated CIR obtained from (2). 
So the ensemble average of hC Rpq li Rpq is 
^ Rp R" R , p , q }=WR hh W" + 

° 2 nW RLS . v W» LSiV (3) 

Where R hh = E{/i p iq h p q }, and 

W = (F p 'F p + <El Lg y 1 F p l F p . 
The PDP of multipath channel within length L g is 
R hh and all the off diagonal elements of R hh are 
zeros . Hence the covariance matrix can be written 
as R hh =diag(p h ). 

where, 

Ph = tPo.Pi p Lch ,0,..,0] r and 

Pi = E{\h pM [n pl l]\\ 
Due to the presence of F p F p in W leads to the 
distortion in R hh . 

So we focus on the method which eliminates the 

spectral leakage of W, instead of calculating W" 1 . 

The covariance matrix of the estimated CIR is 

defined as 

R~ R =WR hh W" 

which can be expressed as follows 

R hK = Z^Wdiag ( Pl Ul )W H , (4) 

Where 

u l is a unit vector whose 1 th entry is one, otherwise 

zeros . Let and t, 

be the L g X 1 vectors defined as 

PK=D g {Rg l )™& 

t, = DgtWdiagiu^W"), respectively, 
where D g (A) is the column vector which contains 
the diagonal elements of A. Then, the relation in 
(4) is simplified as 

PR = Poto + Pih + - + PL g -ih g -i = T Ph 

(5) 

Where T= [ t 0 ,t 1 ,....,t L _J is a distortion matrix 
caused by W. The distortion matrix satisfies the 
condition l[r] ii l>I] 7;ti |[T] iJ | forall i,j, which says 
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that it is a diagonally dominant matrix. The non 
diagonal elements are leakage powers of u ; for all i. 
so that the distortion of W can be eliminated as 

Ph =T l p R = E {g p , q [n p ]}-^w (6) 

PDP at the (p,q) th antennaport on the n p th OFDM 
symbol estimated using the received sample vector 
9 P , q [n p ], 

Where g pM [n p ] = T l D g (R Rp q Kl P , q ) and 
w=T l D g (w RLS , p W R H LSp ). 

PDP ESTIMATION IN PRACTICAL MEMO- 
OFDM SYSTEMS: 

The sample vector received is given as 

9p.qM = D g( h p.d h P.^ +il P.q + e p.d ( 7 ) 

Where 

n pM = T l D g (w RLS , p n q n^W R H LSiP ), and 

e VA = 2Re{T 1 D g (W RLS . p h Ptq n q *W» s , p )} 

assume n pq is an effective noise by AWGN. Then 
the sample average of received sample vector 
BvsM is given by 

^ l r pl p Q 
<5 P .,[np]> W = ^ZZZ 5 ™W 

n p p=lq=l 

(8) 

={D g (h pq hl q )) N + in pq ) N + ie pq ) N , 

Where N is the total number of samples for PDP 

estimation which is given as N = \T p |PQ . 

Where 1^1 = the number of pilot symbols at the 

fe p th subcarrier in a time slot. To obtain the accurate 

PDP N must be sufficiently large. Since 

^D g (h Piq h' p l M )) N ^p h , 

(k p Jn -> °l w and ie p q ) N -> 0. 

It is difficult to obtain such a large number of 
samples at the receiver in a practical MIMO OFDM 
system. With an insufficient number of samples, 
the PDP can be approximated as 

However PDP accuracy with insufficient samples is 
improved by reducing the effective noise as follows 

(g P , q h P ]) N - rtw = (D g (v,/&A + *N 

(9) 

Where, 

z N = (e pq ) N + (n pq ) N - (t£w is a residual noise 
vector, which entry has a zero -mean with N 
samples the error of PDP estimation can be 
calculated as 

e N = ((D g (h p ^ q )) N - p h ) + z N ... (10) 

Since, [p h ] t > 0 for all i, 
P lnlt =^Y.; =1 Y. Q q= ,s Piq [n v \ ,...(11) 



Wheres p q [n p \ is the sample vector of implemented 

PDP estimation with 1 th entry 

s l P , q [n p ] = 

\9p, q [ n p\ ~ l f ' 9p, q [ U p\ > a nW l 
( 0 otherwise 
(12) 

Where 5 <> p ] = [g pA [n p ]] ; and w l = [w] t . 

The implemented technique estimates the average 
of residual noise at the zero-taps of p h , to reduce 
the effect of residual noise z N . The zeo-tap is 
detected at the 1 th entry of p init as 

4 = { 1 if PLu < Ptn ( i3) 

I 0 otherwise 

Where fi th = — 1 p/ nit is threshold value for 

detecting zero-tap. At zero-taps the average of 
residual noise is estimated as 

n Riavg =j-t l 3 = ~ 0 1 Pl nit t l z (14) 

Where 

N z = 51jf 0 t l is the number of detected zero-taps. 
The PDP estimate p h at the 1 th tap afterreducing the 
residual noise, can be expressed as, 

_ \fiinit -KR.avg if Pinit > ^R.avg 

I 0 otherwise 

(15) 

The frequency-domain channel correlation in the 
LMMSE channel estimator is obtained using the 
estimated PDP in (15). 

PERFORMANCE ANALYSIS: 

With imperfect PDP in (15) the LMMSE channel is 
estimated as 

Wf. P = F L D g (p h )F% (F p D g (p ft )F» + aZl kp ) 
(16) 

Where by taking first L g columns of DFT matrix 
F L which is K d X L g matrix is obtained. p h = 
Ph + £pdp is expressed as the estimated PDP. The 
e pdp for the 1 th element is defined as 

E i d = pjvlz- n R>avg if[e N ] t > n R , avg 
pdp ( 0 otherwise 

(17) 

From matrix inversion lemma, F p D g (p h )Fp + 
°n'k ) ' i s converted as 

F p D g (p h )F% + rtl kp y l =A A -A A F p BF H p A" 1 

(18) 

Where A A F p D g (p h )F^ + all kp J l 
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B a Dg (e pdp )(I Lg + F» A-%D g {e pip )y\ 
Then, the coefficient matrix for LMMSE channel 
estimation with p h can be rewritten as 

W f,p = W opt,p + Werr.p. (19) 

Where 

W op t.p = f lD 9 (p h )fp (F p D g (p h )F^ + all kv ) A 
is the coefficient matrix for wiener filtering, and 

W err p = — Fj^Dgip^Fp A' l F p BF p A" 1 
+K D g {.e pdp )F H p F p D g {p h )F H p + a£l) A 

(20) 

With imperfect PDP the error covariance matrix of 
LMMSE channel estimation can be obtained as 

E p = E{ {F L h p q - W f p K LSiPiq )(F L h p q 

-w f ,phs. P . q ) H } 

=(F L - W f p F p ) D g (.p h )(F L - W f p F p ) H 
+a£W f)P F p F v 'W f H pi (21) 

Where K LS p q A diag (X p ) H y p [n p ] 

Using the error covariance matrix the frequency - 

domain MSE of the implemented method is 



MSE f , 



. (22) 



RESULTS AND DISCUSSION: 

The simulation results obtained in MIMO OFDM 
system with two transmit and receive antennas. The 
simulated graph shows the channel estimation by 
taking different parameters. The simulated graphs 
show the MSE performance of LMMSE technique 
using estimated PDP. 



between MSE and channel length. As the channel 
length is increasing the MSE also increases the 
interference get increases. In figure 5, graph shows 




Fig 4: Graph plotted between MSE and channel 
length using implemented MMSE and ID LS 
method 




Fig 3: Graph plotted between MSE and SNR 
using implemented MMSE and ID LS method. 



In figure 3 the graph shows the relation between 
MSE and SNR. As the MSE is decreasing then the 
SNR is increasing i.e., as the mean square error is 
reducing the signal level is increasing which leads 
to increase in SNR ratio. In figure 4, graph plotted 



Fig 5: Graph plotted between MSE and mobile 
equipment speed using implemented MMSE and 
ID LS method 




Fig 6: Graph plotted between MSE and number 
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of samples using implemented MMSE and ID 
LS method 



[9] T. CuiandC. Tellambuia. " Power dela\ pi\> li 1c and noise 
variance estimation for OFDM," IEEE Communication. Lett. 
,vol. 10,Jan.2006. 



how the MSE varies with the mobile equipment 
speed. As the mobile is moving fast then it 
experiences a fast fading in multipath channel 
which increases the MSE. In figure 6 as the number 
of samples increases the MSE decreases. 



[10] H. C. Won andG. H. Im, "Iterative cyclic prefix 

channel estimation for a STBC OFDM system," IEEE Commwi. 
Lett., vol 9, Apr. 2005. 



CONCLUSION: 

The simulated graphs show the comparison 
between the implemented method using LMMSE 
technique and ID least square technique. The 
implemented technique using LMMSE gives less 
mean square error compared to the ID least square 
method. Implemented algorithm holds good, as the 
number of samples increases also observed from 
the result. The implemented technique effectively 
reduces the spectral leakage and residual noise. The 
performance of LMMSE channel estimation 
approaches that of wiener filtering using the 
implemented PDP estimate. 
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